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Abstract 

We consider a family of chiral non-Hermitian Gaussian random matrices in the unitarily invari- 
ant symmetry class. The eigenvalue distribution in this model is expressed in terms of Laguerre 
polynomials in the complex plane. These are orthogonal with respect to a non-Gaussian weight 
including a modified Bessel function of the second kind, and we give an elementary proof for this. 
In the large n limit, the eigenvalue statistics at the spectral edge close to the real axis are described 
by the same family of kernels interpolating between Airy and Poisson that was recently found by 
one of the authors for the elliptic Ginibre ensemble. We conclude that this scaling limit is universal, 
appearing for two different non-Hermitian random matrix ensembles with unitary symmetry. As a 
second result we give an equivalent form for the interpolating Airy kernel in terms of a single real 
integral, similar to representations for the asymptotic kernel in the bulk and at the hard edge of 
the spectrum. This makes its structure as a one-parameter deformation of the Airy kernel more 
transparent. 



1 Introduction 



Random Matrix Theory (RMT) is a rich field of probability theory, with many applications in physics, 
mathematics and beyond. One of the results that exemplifies this success is the appearance of the 
so-called Tracy- Widom (TW) distribution, describing the largest eigenvalue of random matrices, [22J, 
in a large range of contexts, including combinatorics, growth processes and random tilings, as reviewed 
in |23| . The TW distribution function, which can be expressed in terms of a solution to the Painleve 
II equation, is defined as a Fredholm determinant of the Airy kernel. This kernel is universal, in the 
sense that it describes the soft edge eigenvalue scaling limit for a large class of non-Gaussian and 
non-invariant ensembles (see |23j for references). 

In the recent past, non-Hermitian RMT has become an area of intensive research, with applications 
ranging from the fractional quantum Hall effect [9J to quantum chromodynamics [21] . and we refer 
to |18| for a recent review and more applications. The following question thus naturally arises: what 
happens to the largest eigenvalue of a random matrix when it is allowed to move out into the complex 
plane? This question has been answered by one of the authors in [7J for the so-called elliptic Ginibre 
ensemble with unitary symmetry {(3 = 2). This is a Gaussian matrix model generalizing the original 
complex Ginibre ensemble [13] • It can most conveniently be defined as a two-matrix model of a 
Hermitian matrix H and an anti-Hermitian matrix A, depending on a non-Hermiticity parameter 
t G [0, 1), see e.g. [11]. Depending on the way r is scaled when the matrix size n tends to infinity, a 
new interpolating behaviour between the Tracy- Widom (r = 1) and Gumbel distributions was found 
in [TJ. More generally, the limiting two-dimensional eigenvalue point process at the spectral edge 
interpolates between Airy and Poisson, and is described by a new family of kernels, generalizing the 
Airy kernel. 

In this work we address two questions: how universal is this interpolating Airy kernel, and can its 
relation to the real Airy kernel be made more transparent? 

To partly answer the first question we consider the chiral extension of the elliptic Ginibre ensemble 
with unitary symmetry, introduced in |19j . Non-Hermitian models with chiral symmetry can typically 
be formulated as two-matrix models with two complex non-Hermitian matrices (depending again on 
a non-Hermiticity parameter). In the Gaussian case studied here, the model can be considered as 
a non-Hermitian generalization of the Wishart-Laguerre ensemble. In all three symmetry classes 
(real non-symmetric (/? = 1), complex {(3 = 2) and quaternion real (/? = 4)), of chiral non-Hermitian 
RMTs, all eigenvalue correlation functions for finite n can be expressed in terms of a kernel of Laguerre 
polynomials in the complex plane, see [5j [191 0] respectively. We will present a proof that, for j3 = 2 
in the appropriate large n limit, the same interpolating Airy kernel as in [7J is found. In that sense 
the interpolating Airy kernel is universal, appearing for two different symmetry classes of Gaussian 
non-Hermitian models with unitary invariance. This is not unexpected as it parallels the situation for 
Hermitian ensembles, where chiral and non-chiral models with the same /3 share the same soft edge 
spectral behaviour [10] . Our method is similar to that of [7J, and involves the asymptotic analysis of 
a double contour integral representation of the Laguerre polynomials. 

To answer the second question we will show that the interpolating Airy kernel permits a real 
integral representation as an alternative to the complex double contour integral representation given 
in [7J . This explicitly reveals the analogy of its structure as a one-parameter deformation of the Airy 
kernel with the known non-Hermitian generalizations of the sine and Bessel kernels. 

In order to put our results in a more general context, let us briefly review the different scaling 
regimes in non-Hermitian RMT, and what is known about their universality. The first region that was 
studied was the so-called bulk of the support of the limiting eigenvalue distribution. Here different 
behaviours occur depending on how the non-Hermiticity parameter scales with n; in the "weakly 
non-Hermitian" regime, n(l — r) =: a 2 — > a 2 G (0, oo), a transition takes place between "strongly 
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non-Hermitian" statistics (a^ — > oo) and what we will call "essentially Hermitian" statistics [a\ — > 0), 
described by the sine kernel, see [1 1|, 112] . Typically, in the strong non-Hermiticity limit, the correlations 
can be obtained from maximal non-Hermiticity (r = 0) by a simple rescaling of the eigenvalues, see 
e.g. 1 121. In the strong non-Hermiticity limit, universality in the bulk was proved rigorously for 
general normal random matrices in the (3 = 2 case in [8J. It remains an open problem how to amend 
the Riemann-Hilbert approach, which has proven so successful in Hermitian RMT, to models with 
complex eigenvalues, see [16] for a discussion on this. At weak non-Hermiticity a heuristic argument 
in favour of universality was given in [12] and [2] for iid matrix elements and weight functions with 
non-Gaussian harmonic potentials respectively. 

Chiral models show a different (hard edge) spectral behaviour at the origin. It depends on the 
finite multiplicity v of the eigenvalue at the origin, which is kept fixed when n goes to infinity. Again 
the same transitional regime of weak non-Hermiticity occurs, as identified in [3], following [llj. At 
the origin in this regime, universality is known to hold in the sense that the same Bessel density in 
the complex plane is obtained from two different Gaussian one- and two matrix models, see [20J and 
[19], respectively. 

At the edge of the spectrum close to the real axis which is our concern here, the appropriate 
weakly non-Hermitian scaling limit was identified in [7] (see also [33]). Here the transition occurs 
in a different regime, o\ = (1 — r n )n 1//3 — > (0,oo), when the bulk and hard edge statistics are still 



strongly non-Hermitian. Our work is the first argument presented to suggest the universality of the 
corresponding interpolating Airy kernel, as we obtain it from a second Gaussian model in a different 
symmetry class. 

Let us also mention the decay of the eigenvalue density at the spectral edge in the limit of strong 
non-Hermiticity. For all three symmetry classes of the elliptic Ginibre ensembles (/? = 1,2,4) the 
density decays as a complementary error function, as was shown in [18] (and references therein). In 
that sense this behaviour is also universal. Note that this is consistent with the Poisson statistics of 
the extreme eigenvalues, which is found on the microscopic scale and shifted out from the spectral 
edge, as opposed to the mesoscopic scale where the decay of the eigenvalue density appears. 

Our paper is organized as follows. In Section [21 we briefly introduce the model and summarize our 
results, which are: An elementary proof for the orthogonality of the Laguerre polynomials on C which is 
detailed in Section [3] (for the corresponding result for Hermite polynomials, see OHE]). An alternative 
representation of the kernel of Laguerre polynomials as a double contour integral shown in Section [H 
which serves as a Lemma to prove the convergence in the scaling limit of the eigenvalue process to the 
interpolating Airy process in Section Finally, we provide a simple real integral representation of the 
interpolating Airy kernel as an alternative to the previously known double contour integral form; the 
proof is given in Section [6J Here we also check that the complementary error function density decay is 
obtained from the interpolating Airy kernel in the appropriate strongly non-Hermitian scaling limit. 

2 Summary of results 

We consider a chiral two-matrix model in the unitary symmetry class which is the non-Hermitian 
extension of the chiral Gaussian unitary ensemble. For any non-negative integer u, let P and Q be 
n x (n + u) matrices with iid centered complex Gaussian entries of variance l/(4n), 



where dP denotes Lebesgue measure on the space of complex n x (n + v) matrices (identified with 
^2n(n+u)y -^Yg are i n t eres t ec j i n th e eigenvalue distribution of the (2n + u) x (2n + v) random Dirac 




dIV(P, Q) = a,n exp -2nTr(PPt + QQ ] ) dPdQ 



3 



V 



(2.1) 



matrix 
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where r £ [0, 1] is a non-Hermiticity parameter. Equivalently, for the distribution of T> is given 

by the probability measure 
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(2.2) 



with d<£ and d^ denoting Lebesgue measure on the spaces of complex matrices of appropriate dimen- 
sions. The spectrum of T> consists of an eigenvalue of multiplicity v at the origin, and 2n complex 
eigenvalues {±z/ c }^ =1 which come in pairs with opposite sign. Without loss of generality, we choose 
the Zk to have non-negative real part. With probability 1, the z k are all distinct and have strictly 
positive real part, and we consider this situation from now on. Changing variables and integrating 
out the eigenvectors induces a probability distribution on {z k }^ =1 , 



d¥^( Zl ,...,z n ) 



where d 2 z = dxdy denotes Lebesgue measure in the plane and Z 
we have introduced the Vandermonde determinant 



An(^ 2 )| 2 [J < U (4)X{R CZk >o}d 2 z k , (2.3) 

is a normalizing constant. Here 



k=i 



n 



A n (z) : 



l<j<k<n 



and the real and positive weight function 

<'"(*) ■ 



zr +1 exp 



2rnRe(z) 



(1 



2n\z\ 



(1 



(2.4) 



where K v is the modified Bessel function of the second kind (also called the MacDonald function). 
The limiting case r = 1 reduces the two-matrix model to the standard Hermitian chiral Gaussian uni- 
tary ensemble, which is equivalent to the unitary Wishart-Laguerre ensemble of non-negative definite 
Hermitian matrices PP^ with non-zero eigenvalues {z|}" =1 . For details on the model, we refer to |19j 

with n = v / (l-r)/ v / (l + r), V H = W, a = 2/(1 +t),A = VTT^P, B = iy/TT^Q, C = t$ and 
D = M. 

The eigenvalue measure (12. 3p has the structure of a determinantal point process and as was shown 
in |19j . the model can be solved by introducing Laguerre polynomials orthogonal in the complex plane 
with respect to the weight function (I2.4I) . The orthogonality of the Laguerre polynomials on C was 
conjectured in [3], however, the proposed weight is only correct for v = 1/2. Let L v - denote the 
ordinary generalized Laguerre polynomial with parameter v of degree j, defined by the orthogonality 
relation 

— o jk . 



L u Ax)L v k {x)x v e- x dx 
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Proposition 1. For any integers j,k,u > and real numbers a > b > the following orthogonality 
relation holds. 

r L» (cz) L u k {cz) \z\ v e bKe ^K u (a\z\) d 2 z = h%6 jk , (2.5) 



MM 



with squared norm 



andc: = (a 2 - b 2 )/{2b). 



Tjjj + v) 
a jl 



a\2j 
b 



2a 



b 2 



v+1 



(2.6) 
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In stating Proposition Q] we have changed back to unsquared variables over the full complex plane 
and general real parameters a > b > 0. The case of interest to us, (|2.4p . follows by setting a = 
2n/(l — r 2 ), b = 2rn/(l — r 2 ) leading to c = njr. While (|2.5p and f)2.6f) were given in [19] as "easy 
(to) verify", no formal proof was offered. A first proof of (|2.5p was given in [4], Appendix A, by 
explicit computation of the moments f d 2 z Wn U \z 2 )L v - (cz 2 ) z 2k in terms of hypergeometric functions. 
A more elementary proof of Proposition [T] is presented in Section 

An alternative proof could be given exploiting the fact that in this symmetry class the orthogonal 
polynomials P n {z) are directly given by the expectation value of a characteristic polynomial P n (z) = 
E[det(X> — zl)], see [6]. The latter can then be computed explicitly for the Gaussian two-matrix model 
(12. 2p using Grassmannians, along the lines of [5]. 

In the limit b — > (r — > 0) the weight w(z) in (|2.4p becomes rotationally invariant and the proof 
of Proposition [T] is straightforward, replacing the Laguerre polynomials by their leading monic power 
(and multiplying (|2.5p with the appropriate power b 2 i). 

By a standard argument, performing suitable row and column operations on the Vandermonde 
determinant, all fc-point correlation functions of the point process (|2,3p can be expressed as deter- 
minants of a correlation kernel defined in terms of orthogonal polynomials with respect to Wn' u (z)/\z\, 



(^-n' T (Ci> Q))l<i,j<k 



,eCj>0} 

3=1 



■ det 



where, by (|2.5p . 

^(d,C 2 ) = K-^K^Cl)) 172 ^ (0y L « (^) L l ) ( 2 - 7 ) 

if Re(j > 0, j = 1,2, and K^n (d, C2) = otherwise. When taking the large n scaling limit at the 
edge of the spectrum the following representation will prove very useful. 

Lemma 2. For any integers n > 1, a > 0, parameter r £ (0, 1) and complex numbers z\ 7^ and Z2 
t/ie identity 



— ■ r 2fc /t! ra/ „ , rQ/ „ , _ r 2n e 22 / J / J ^(u - l)^ a e 21 !^ 1 ) 22 (^T 

«!!/ (r 2 u — u)(i> — l)(u — 1) 

(2-8) 

holds. Here 71 is a simple closed contour encircling the point u = 6wt not n = 1, w/uZe 72 encircles 
both the point v = 1 and i/ie entire contour 71 in sitc/i a way £/ia£ r 2 w - for u 6 71 and v 6 72 . 

Now we state our main result, giving scaling limits for the eigenvalue distribution of T>. We 
emphasize that the following result holds for any fixed v. In principle, our method could be extended 
to cover the case when v scales with n (see Remark [6]) , the difficulty here is to identify the global 
spectral distribution, a non-Hermitian generalization of the Marcenko-Pastur law. 

Theorem 3. Let v be a fixed non-negative integer. Given a sequence {r n }^ =1 C [0,1), put a n = 
^n) 1 / 6 -^/ (1 — T n ). Let Zj = Xj +iyj, j = l,...,n be the eigenvalues with positive real part of 
the random Dirac matrix T>, forming a determinantal process with kernel KX% Tn ■ For the choices 
of scaling parameters a n , b n and c n specified below, put Xj = (xj — c n )/a n and i/j = yj/b n , de- 
fine the rescaled two-dimensional eigenvalue point process Zn n ' u = {(xj,yj)}^ =l , and let Fn"' u (t) = 
¥^" u [maxi<j< n {ij} < t] be the last particle distribution of Zn"' u - 
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(i) Suppose o n — ^ o 6 [0,oo) as n — > oo. Choose a n = (2n)~ 2 / 3 , 6 n = cr n (2n)~ 2 / 3 and c n = (l + r n ). 
T/ien Z^"'^ converges weakly to Z^, a determinantal point process on C = R 2 u>i£/i correlation 
kernel 

Iff e -^u+r) 1 ) 2 + ^u 3 +i^ 1 u-\(av-r) 2 ) 2 + iv' i +ii2V 



C(Ci,C 2 ) : = -r^ / / dt; — — , (2.9) 

47r 5 /2 J R+ . S J R+ . S l(u + V) 

where Q = £j + irjj, and Fn n,u (t) converges to the last particle distribution F a (t) of Z^. 
(ii) Suppose a n — > oo as n — > oo. Put f n = (1 + r n )/2 and choose 

:(2n)- 2 / 3 , (2.10) 



and 



V6 log a n 

5/2 

bn = f -V4 g» , (2n)- 2 / 3 , (2.11) 
c n = (1 + r n ) +a„ (31oga n - ^ log(61oga n ) - log(2vrf 3 / 4 ) j . (2.12) 



Then Zn"' u converges weakly to a Poisson process on M? with intensity 7T 1 / 2 e ^ v , and Fn n ' v (t) — 
e~ e , the Gumbel distribution. 

The interpolating Airy kernel (j2.9j) was derived as the scaling limit for the corresponding non-chiral 
Gaussian ensemble in [7J. We show that this kernel has an alternative real integral representation. To 
clarify the analogy with the transitional Bessel and sine kernels below, we consider the scaling limit 
corresponding to b n = a n = (2n) _2//3 rather than b n = a n a n = <r n (2n) -2 / 3 in the first part of Theorem 
[3j for a > 0, define the rescaled version 

C2) + 6 + 

= ^ / d« / dw (2.13) 

of the kernel (|2Uj) . 

Proposition 4. T/ie interpolating Airy kernel has the real integral representation 

p -^?+^ 2 )+^ 6 + #(Ci+C 2 ) roo 

C2) = ^ / dt e tCT Ai (Ci + a 4 /4 + t) Ai (C 2 + a 4 /4 + t) , (2.14) 

where Q = £j + irjj and Ai is f/ie Airy function. 

The convergence of the integral in (|2.14p can easily be seen from the large argument asymptotics 
of the Airy function, 

Ai(*) = ^*^(l + (1)), (2.15) 

as z — t- oo with Arg(z) bounded away from n. 

In the limit a — > the integral in (|2.14p becomes elementary (see (2.37) in [T7j for this integral 
form) , 

^(Ci,C 2 ) ^— / dt Ai(fr + *)Ai(& + t) = ^/c A (ei,6), 

V 71 " Jo V 71 " 
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where K. A is the Airy kernel. Proposition |4] completes the picture of the known non-Hermitian gener- 
alizations 

i2-( 7 7i+ 7 7i) /•! 

JCf (Ci,C 2 ) = 6 ^ 3/2 / dt e"^ cos(t(Ci - C2)), (2.16) 



121 of the sine kernel )C S and 



[3lll9| of the Bessel kernel JC B ' U , in the bulk and at the hard edge of the spectrum, respectively. These 
are obtained in the limit when a = lim n _ 5 . 00 \J (\ — T n )n £ (0, oo), and are both given as simple one- 
parameter deformations of the ordinary sine and Bessel kernels, which are recovered in the essentially 
Hermitian limit, 

~rS( t t \ *=$ e^^+^)sin(vr(6-6)) e~W,+nl) 

^(6,6), 



7T 



where we have rescaled the imaginary part with a factor a to get a finite limit, analogously with 
the definition (12. 9ft of K.^. Note that the bulk and hard edge transitions occur in a different regime 
(1 — r n ) >c re -1 compared to the soft edge (1 — r n ) X n -1 / 3 . 

Remark 5. From the representation (]2.14p it is also transparent that our results are consistent with 
the expected error function decay of the eigenvalue density at the spectral edge when n — > 00; the 
density of rescaled eigenvalues 2j := y/2n(l — T n )~ l / 2 (zj — (1 + r n )) is given by 



»»«> : = (P + + P + + J - 2^(1+7) e *<«- 

(2.18) 

where £ = £ + i-Vi T = um n-»oo r «> an d erfc(n) = dt e * is the complementary error function. 

See Section [6] for a sketch of this argument. 

3 Proof of Proposition [I] 

In this section we present an elementary proof of the orthogonality of the Laguerre polynomials in 
the complex plane, Proposition [H using induction in v. We will make use of known recursions for the 
Laguerre polynomials and modified Bessel functions which are of depth 2, so we first show the base 
cases v = and v = 1. 

To that end we will use an integral representation for the modified Bessel function of the second 
kind, ([I5|, 8.432.6), 



™-i®T^ m - (3 - 1) 



as well as a complex contour integral representation for the generalized Laguerre polynomials 

1 f e -xu/{l-u) 

L v Ax) = — 6 du- . (3.2) 
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where the contour 7 encircles the point u = but not u = 1. We will choose 7 to be a circle 
of radius r < b/a, centered at the origin. For v = 0, (|3.2p is given for instance in |24| . equation 
13.53. Alternatively this relation can easily be derived using the generating function for Laguerre 
polynomials. The case v = 1 (and indeed any positive integer u, although not needed here) then 
follows easily using the identity L^ +1 (x) = —L^ +1 (x). 

Inserting the representations (|3.ip and (|3.2p with v = into (|2.5p gives 

<L° Ll) = f dx r Ay e- f *W * / d,^f / tte? -"- , »»»'"-»» 



J ' ' /-co J-00 Jo 2t (2m) 2 /-, J, (1 



1 ^.00 /-oo /'oc y y e _iC-i r (a ; +tA) 2 -f 7 (y+t J B) 

dx dy dt 



4>&u$dv m ^-^ ^ -fc+1 , (3.3) 



where 



2(27ri) 2 y_ 00 y.oo y ; A A 



^4 = 4 I c( — + 



I? 



a 2 ^ ^ 1 — w 1 — ?j 
2ic / u u 



a 2 V 1 — ti 1 — f 



and 

' _ 4 v ^ ' " ; ~ a 2 (l-u){l-v) \* t- 2 
depend only on the variables u and v. Note that Re(C) > for any -u, v on 7 since |u| = \v\ = r < 
6/a < 1, by our choice of 7. This means that we can apply Fubini's theorem to calculate first the 
Gaussian integrals in x and y and then the convergent t-integral. Inserting all this into (|3.3p gives 

2vr 1 J" du f dv 1 
27r / a 2 \ fe 1 /" du 



/rO r0\ 27F 1 I dl1 I 



(a 2 - 6 2 ) V 6 2 / 2m J 7 

n / 2 \ k 



(a 2 -6 2 ) V&V ^ 



-hjSjk, 



proving the base case v = 0. 

The result for u = 1 follows along the same lines and we only give some intermediate steps. A 
straightforward calculation using integration by parts shows that, as long as Re(C) > 0, 

/-> />f ^- c '- ^ " w (? + ^) " <-> 

implying that 

„, ,,, a 1 I™ 1 I™ , f°°i Ii Ii (i'+S ! )e" IC "I |,t ""'-i (,+ ' B)! 

^- £ »> = 4(2^ L ix iJ v L it t iu t iv W^W^TW^ 

4an 1 f du f dv 1 

= (a 2 - 6 2 ) 2 (2vri) 2 % ^+1 J 7 ^+1 (1 - {a/b) 2 uvf 



4a7r 



(a 2 - 5 2 ) 2 
hjSjk, 



2 \ 3 
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since 



dv 



2m J v k+l (l - {a/bfuvf 



Res„=o v~ 



OO / 9 \ 7 

a u N 



— ' aru^ 



(-w) = (* +1 >(-*r 



9 \ ft 



6 2 



Next we can make the induction step using the recurrence relation for the modified Bessel functions 
of the second kind, ([15]. 8.486.10), 



xK v+1 (x) = xK u ^i(x) + 2vK v (x) , 



(3.5) 



as well as the following two recurrences for Laguerre polynomials, ([IS], 8.974.3, 8.971.2 and 8.971.5 
respectively) 



Ln +1 ( x ) — ^ ^™ 



m=0 



and 



(3.6) 



(3.7) 



m=0 



Suppose (I2.5|) holds for a fixed v and v — 1 and every j, A;. Without loss of generality we consider 
j > k. By equations (|3.5j) through (j3.7|) and the induction hypothesis, 



: y d 2 z e h ^L) +1 ( cz ) (cz) f ^i/^^ (o|z|) + 
Ic 



z\ 2 \z\ v - 1 K v ^ (a\z\ 



d 2 z e bRe ^ 



-[Y^LUcz)) (J2 L n(cz) )\z\ v K u (a\z\ 



+4 Ut^H-O' + l^lN (^^(^-(fe + l)^;}^) ) ^r- 1 ^-! (a\z\ 



\ m=0 

j k 



n=0 



= V E E^™' L ^ + 72 E E^™ S L n X ) " E^J+l ' ^ 



m=0 n=0 



m=0 n=0 



n=0 



/ r 1 r f— 1\ i / T U—l T V— 1\ 



m=0 



2;/ 



n=0 n=0 



(3.8) 



Note that the last two terms contribute only when j > k and j = k respectively. Inserting the 
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expressions (12, 6p we obtain for the two sums 

2f\^ f „ v 2 ^ iv i 2i/ A vr(n + ^)! a 2n / 2a 



a ' — ' '" c ' — ' " a < — ' a n! 6 2n \ a 2 — b 2 

n=0 n=0 n=0 



zv 2 4o 2 A 7r(n + i/- 1)! a 2 " / 2a 
+ (a 2 - 6 2 ) 2 ^ a~^! "o 2 ^ la 2 -6 2 

v ; n=0 v 

47ri/(2a)*' A (n + v-l)\a 2n .. . 2 , 2 . 
a(a 2 - 2 )^ 2 E n! + ^ ~ nh ) 

V ' 71=0 

A-Kv{2a) u {k + v)\a 2k+1 
{a 2 -b 2 ) u+2 k\ W 
(fe+j> x 

c 2 n k+l 



after telescoping the sum. Thus (|3.8p becomes 

, > 2 

r- 



itv+1 rf+l\ _ (^+ 1 ) t/ ^-l/i c-i j; \ ^ i + 1 ) 2 ^-l^ 
( L j > L fc /- ^ "fc+iU - i 1 ~ d jU) + 72 ft fc+l d J* 



i-^^ + i + ^Oi* 

7r(fc + t/+l)! (a\ 2k ( 2a \ u+2 x 
a hi Kb) \a 2 -b 2 ) jk 



which concludes the proof. 



4 Proof of Lemma [2] 

In this section we prove the alternative complex contour integral representation for the kernel, Lemma 
[21 To that end we will derive two different contour integral representations for the Laguerre polyno- 
mials, and insert them into (|2.8p . This will obviously break the symmetry in the arguments z\ and 

The first integral representation follows starting from [I], 22.10.7, 

L%(z) = -— <p dw-^ ^j-r. (4.1) 

valid for any integers k, a and z ^ (for simplicity we exclude non-integer indices a to avoid branch 
cuts). The contour 7' includes the pole at w = z but not the origin w = 0. 

Now make the change of variables v = w/(w — z), dv = —z/(w — z) 2 dw, giving rise to 

p z r „.k+a p —zv/(v—T) 
L?(z) = — & dv , , . (4.2) 

Together with the sign from the Jacobian and the inversion from the Mobius transformation, the 
contour of integration is mapped to a positively oriented contour encircling the pole at v = 1 but not 
the origin v = 0. Without loss of generality it can be chosen subject to the second condition in Lemma 
[2]and thus taken as 72. 
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To prove the second integral representation, 



\a-l 



^- M/M > (4-3) 



(fc + a)! 2vri(-z) Q J 71 it 
with 71 as in the statement of the Lemma, we use induction over a. The base a = is given by (|3.2 

1 C p -zu/{l-u) 



2-rri J lt (1 — u)u k+1 
Using the known recurrence relation for Laguerre polynomials (e.g. [15], 8.971.4) 

zLl + \z) = (k + a + l)L%(z) - (k + l)L% +1 (z), (4.4) 

we can establish the induction step. Suppose (14. 3D holds for a fixed a and every k. Then for any k > 
and z ^ 

k\Ll + \z) k\L%{z) (k + l)\L<* +l (z) 



(k + a + l)\ z(k + a)\ z{k + a + l)\ 



1 / J »"-,f'e-/'-)f 1 -i 



2mz(-z) a f 7l u k+1+a 

1 I &u { }~ Ur , e- zul{1 - u) 



2<Ki(-z) a + l J 7l u fc+i+«+i 

where we used the induction hypothesis in the first line. 
Combining the two representations (|4.2p and (j4.3f) gives 



/ L „ / ' ^("- jo y y- 1 1 (^v\ k e- z2 —- zi — 

I \ (q, — 1 \q, I „ \ „ 



(2vri) 2 zf J 71 J 72 \{v-l)u) f^ Q u\uJ (v-l)(l-«) 

(4.5) 

After computing the geometric sum, 

El / r^vA 1 / t z v \ 1 . 
- = 7-5 r - r-o r, 4.6 

u \ U J [T Z V — U) \ U J [T Z V — U) 



we only have to show that the second term in (|4.6p does not contribute to (|4.5p in order to arrive at 
the desired result (|2.8p . To see this, we change variables v = y/(y — 1), di> = l/(y — l) 2 dy, thereby 
mapping 72 to a new contour -y' 2 which will not enclose the pole at y = u/(u — r 2 ) since 72 does 
by definition, and the Mobius transformation v 1— >• y involves an inversion. The second term of the 
integral in (|4.5p with respect to y then becomes 



e 22 / y a e~ Z2V 

dv [ , " 7T ) 7^5 T7 -r = <p dy 



(v-l)J (t 2 v-u)(v-1) L, y {T^y/(y-l)-u)(y/(y-l)-l)(l-yf 



dy , » - r- r =0 



£ 



y<* e -Z2V 



<i 2 {y{T 2 -u)+u) 
by Cauchy's theorem. 
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5 Proof of Theorem [3] 



To prove our main result, we show that (a kernel equivalent to) the correlation kernel fC% Tn of the 
eigenvalue point process, appropriately rescaled, converges point-wise to the limiting kernel associated 
with the sequence {r n } in the statement of the theorem. The convergence of processes and last particle 
distributions then follows from a dominated convergence argument, see [Tj- We will assume that r n 
is bounded away from 0; the case r n — > reduces to the (simpler) radially symmetric case where the 
kernel can be expressed in terms of monomials rather than Laguerre polynomials; our results remain 
valid in this case. The analysis is very similar to what was done in great detail in [7] for the case 
of the elliptic Ginibre ensemble, so for clarity we will skip some technical details and refer to those 
calculations to give an idea of how the present analysis can easily be made completely rigorous. Our 
strategy is to use Lemma [2] to represent KXf and find its large n asymptotics by steepest descent 
analysis of the integrals. 

It is known that the mean eigenvalue density of T> is asymptotically constant in the ellipse 
{(£.iV) '■ £ 2 /(l + T n) 2 + f? 2 /(l — T n ) 2 < l}. We are interested in the scaling limit around the spec- 
tral edge on the positive real axis, so we define the new scaling parameter S n := c„ — (1 + T n ) 2 <C 1. 
The choice of parameters will be determined by the analysis, however in order to estimate orders of 
magnitudes of error terms, we will assume when needed that the parameters scale correctly with n, 
that is 



^_„-2/3 



and 



and 



VTog CFn 

5/2 

(loga n )V4 



5 n x a n y/\oga n n 2/3 

if u n — )• oo and a n x n~ 2//3 , b n >c a n n~ 2 / 3 , 5 n = otherwise. 

Using Lemma [2] to represent the kernel (|2.7p and the large argument asymptotics of the modified 
Bessel function of the second kind ([I], 9.7.2), 

K v {z) = ^f z e- z (l + o(l)), 

the kernel Ain' Tn ' (Ci > C2) := a n b n ICn T " (zi, Z2) of the rescaled eigenvalue process Zn n ' u can be written 

y / \zi\\z 2 \n 3 / 2 a n b n f T n \z 1 \\z 2 \ 



TT 



5/2 



cxp 



/ z 2 (\z 2 \ + \z 2 \-r n Re(z 2 + z 2 )) 
n 12 log r n + ^ 

I T n (1 - T 2 ) 



dw l( p dw 2 21 1 i ) — 2 , (l + o(l)), (5.1) 



'7i J 72 

where Zj := c n + a n £j + ib n r]j, 

f atj := Re(z 2 ) - c 2 = 2c n a n ^ + a 2 n i) - b\j] 2 
\ (3j := lm(zj) = 2b n rij(c n + a n £j), 



(5.2) 



and where we have defined 

fnH = ---^ + ^gw, (5.3) 

T n W-l 
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with the branch cut of the logarithm chosen along the positive real axis, in order that f n be analytic 
in neighbourhoods of its saddle points. 

The dominating contribution to the integral will come from saddle points of /„, given by f' n {w) = 0, 

or 

w 2 +( C "~ 2 M w + 1 = 0, 



which has two real solutions, W- < — 1 < w+ < 0. 

Remark 6. Note that, more generally, if v scales with n, v = an, the integral can be evaluated in 
terms of saddle points of the function g n (w) = f n (w) + a(logii; — log(w; — 1)). In principle the analysis 
can be carried out for this case as well, however the edge of the spectrum c n then moves on a global 
scale and is only implicitly determined, which makes the calculations more difficult. However, the 
present analysis does apply for any fixed v. 

Explicitly, the two saddle points are given by 



or, when expanded in 5 n /(l — r n ) 2 <C 1, 

= ~ ~ _ f 1 " rl(l + r n )—^-,(l + o(l)) j (5. ! ) 



T n T n (l - -J" 2 ) V "" (1 -T n , 



and 



respectively. Since 



'"+ = - T " + TT^i) (' - (1 + Tj O^' 1 + o(1)) ) (5 ' 5) 

2cl 1 

/nW = 7 ^TT3 2 ' ^ 5 - 6 ) 

T n (w — iy w z 



we find that 



and 



/•//, , , r n(l T n) ( \ 2(1 T n T%) 5 n \ 



/ " < - +) - (' + 2(1 (i + :;:)^(i-\)' (i + < °< <5 - 8) 

whereas \ fn(w±)\ X 1. Note that, here and in the following, correction terms depending on 5 n are 
identically zero unless a n — > oo. 

We will choose the contours so that 71 passes through w + and 72 through W- and are both 
tangent to vertical lines, giving the directions of steepest descent, by (|5.7p and (|5.8p . Computing the 
asymptotics of the kernel, the main contribution will come from neighbourhoods of the saddle points, 
so neglecting small contributions away from the saddle points we can take the contours of integration 
as the lines w + + iR and W- + iM. respectively and change variables w\ — w + = ix, W2 — W- = iy. By 
expanding in powers of x and y, we can find the asymptotics of the integral in (15. ip . Before proceeding 
with the saddle point approximation, we calculate asymptotics of the constant and subleading terms 
at the saddle points. 

We note immediately that, provided r n is bounded away from 0, the factor 

w 2 (wi - I) 



(W2 - 1)U>1 
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in the integrand is approximately constant, taking the value 

, V u - =T n "(l + o(l) 



which cancels the pre-factor 



2 



so that the limiting kernel will be independent of v. Furthermore, 

(wi - 1)(W2 - 1) = (1 + Tn)2 (l + o(l)), 
and the remaining factor in the denominator becomes 

r n w 2 -w 1 =T n x-y- — (1 + o(l) . 

(1 - ^n) 

The leading order term of the exponent is 

/„(»-) - /„(»+) = -2.ogr„ - 2„ - fjl^ - (1+T J (1 _ Tn) a + 

We will also need expansions of the sub-leading terms, 

w W- 1 „, 9 . 



2 



where 



and similarly 



(w - 1) (to- - 1) - 1) 



+ n i _ %a7i — rri 1 + °( 1 ))' 



(«;_-l) l+r n (l + r n ) 3 (l-r, 



■(l + o(l)). 



(W+-1) l+r n (l+r n ) 3 (l-r„ 

Collecting all the terms independent of the variables of integration from these calculations, and 
the exponential pre-factor, gives 

z\ \zi\ 2 + \z 2 \ 2 -T n (Razl + Razl) w + (ax + iPi) W-(a 2 -if3 2 ) 

2l0gTn + ~ n (T^f) + ^ - UW+) + r n (w + - 1) " r B («—l) 

Mai + ° 2) -d + o(D) - (1 + OK)) - ,(1 + o(l)) 



(l + r n )3(l-r n ) v v " 2c*(l-7*) v (l + r n )3(l- Tn ) 

+ Cn(/3i-/3 2 )-2vri + o(n- 1 ), (5.9) 

for some constant C n . The correlation kernel -Mn' Tn (Ci;C2) can be replaced by any kernel of the 
form e 9n ^~ 9 "^A4n Tn (Ci> C2) without affecting the underlying point process since the determi- 
nants giving the correlation functions remain the same; we will thus consider the equivalent kernel 
e~ Cn ^ /3l ~ l3 ^Mn Tn ((i,(2), which we still call Mn Tn with a slight abuse of notation. Inserting (j5.9|) into 
(|5.ip gives 
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•Mn ,T "(Cl,C 2 ) 
T n n 3 / 2 a n b n e (i+-n) 3 (i-r„) 

= 7T 5 /2(l +Tn )3/2y(r3^J 



cxp 



-n 



5 n («i + a 2 ) 



+ 



(^i 2 + /3 2 2 ) 



(l + r B ) 8 (l-r„) 24(1 -r 2 ) 



dx / dy- 



B „(/ n ( w _^-/ n ( w _)- /n ( w+ ^) + / n(w+ ) + ^i^-^^) 



where 



and 



i(s - r 2 y) + 

2<5 n 



{w^-iy (l + r n ) 2 V (l + r n ) 2 (l-r n ) 



1 



1 + 



2r /) 
^ 'n u n 



(l + r n ) 2 (l-r n ) 



■(l + o(l)), (5.10) 

(5.11) 
(5.12) 



(W + - l) 2 

5.1 Case a n — » a < oo 

In this case we may choose 5 n = 0. This means that, in order to get a convergent integral in (|5.10p . 
we have to shift the contours of integration slightly into the complex plane. For the details of these 
considerations, we refer to the very similar arguments in [7]; here we will formally integrate over the 
real line to clarify the presentation. Since r n tends to 1 in this regime, 

f:(w±) = + (i)) = —a + (i)). 

After a change or variables, y = 2(2ra) _1//3 t, x = 2(2n) _1//3 s, in the integral in (|5.10p . (|5.7p and (|5.8p 
give 

n(/ n (™_ + iy) - f n (w-)) = -\a 2 n t 2 + U 3 + 0{ n - x lH% 

and 

n(f n (w+ + ix) - /„(«,+)) = -\a 2 n s 2 - % -s 3 + ©(n" 1 / 3 /). 
In the remaining expressions, we may directly insert r n = 1 and (I5.10P becomes 



^' T "(Cl,C2 



exp 



a n 6 n n 4 / 3 



2 2 / 3 vr 5 /2 a , 

ds / dt 



n 



16(1 - r„) 

e -^2(t 2 +S 2 )+f(t 3 - S 3 )+i2- 4 /3 n 2/3 ((a2 _ i/32)t _ (ai+i/3l)s) 

i(s - i) 



l + o(l)). 



Recalling the definition (|5.2p of our auxiliary variables Oj and and choosing a n = (2n) 2 ' 3 and 
6 n = a n (2n)" 2 / 3 

p -|(cr„t-r ? 2) 2 + |t 3 +i6t-5(o-nS-»7l) 2 -|s 3 -^l« 

(i + (i))->JC^(Ci,C2) 



ds / dt- 



gives 

-^n' T ' l (Cl,C2) = -47 2 - , ' / - , 

47T 5 /2 J R J R l( S - t) 

as claimed. 

For the weak convergence of point processes and convergence of the last particle distribution to 
follow, it now suffices that there for every £o is a uniform bound B a £ (Q > \Mn T "(CiC)\ whenever 
Re(C) > Co ; such that Jjw^)-^ d- 2 CB a £ (() < oo. It can be shown that there is such a bound of the 

form B a £ (Q = C a £ e~2^ +V \ The argument parallels that of [7J, and the details are omitted here. 
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5.2 Case a n — > oo 

This case requires a subtle choice of scaling parameters, to be determined by the analysis. Change 
variables, t = nf%(w_) y and s = y/—nf^(w + ) x. The imaginary part of the denominator in (|5.10|) 
becomes 

x - r 2 n y = iTn J^^. {l + (l))(s - t), 

so using the asymptotics (15, lip and (|5.12j) . and noting that 

1 1 / 1 \ 

i-W-; — v? (i + o(i)), 



24(1 - rl) 2(1 + r n )3(l -r n )\ n (1 + r n f 
we can write (|5. 10[) as 

AC T "(Ci,C 2 ) 



na n b n e (i+^a-^n) 

- exp 



7T 5 /2(l_ r , 



2> 



-n| ^ + ^) + (fl+g) u _. 



(l + r n ) 3 (l-r n ) 2(1 + T re ) 3 (l - T n ) V (1 + T, 



where 



An := — / f „, si T^^ 1 - iai )> 



/ Ti 

B n := ==% 7-77(^2 + *a 2 ), 



r nV ^I)(u;_ - l)^ 



and 

Pn ■= , 2SnVE >0. 

We note that p n x \/log <r n 3> 1, by assumption on cr„. Expanding, we find that 



/77- 

Im (A n - fl n ) = -(ai + a 2 ) * — (1 + o(l)) « p n , 

V(l -T„)(l +T„)'V 2 

and 

Re (4„ - B n ) = 00! - ft) ^ — (1 + o(l)) x af^log^) 3 / 4 » p n , 

for /3i ^ /?2- To calculate the integral in (|5. 13|) . we write the denominator as an integral and then 
compute the Gaussian integrals in s and t, giving 

ds / d t — 7 — = da dt du e -\^ + s 2 ) + A nS+ B n t-u {Pn M{ S -t)) 

p n + z(s - i) J M J R Jo 

= 2ir du e \( A ^~^) 2 +\(B n +iuf-u Pn 
Jo 

poo 

= 27re^ (j4 " +B " ) / du e -^ 2 -^(Pn+i(A n -B n )) 



l + O li)) (5.14) 



(p n +«'(A» - #n)) V VPn 
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where the last integral is calculated using integration by parts and a suitable splitting of the integral. 
Inserting ()5.14j) into (|5. 13|) leads to some cancellation, and we are left with 



AC' T ' l (Ci,C 2 ) 





y/(l + T n )^Ea n b n e (i+^(i-r n ) 

: 7r 3/2 v ^r^rj (5n( i + _i_( A _ A)) 



x exp 



-n 



5 n (ai + a 2 ) 



+ 



(l + T 2 )6 n (P 2 + P 2 ] 



(5.15) 



.(l + r n )3(i_ Tn ) 2 (l + r n ) B (l-r n )3 

where we have conjugated away a factor of the form e^™^ 2 ^ 2 "" 1 ^ 1 ) coming from (|5.14p . In terms of 
the original variables, the exponential terms of (|5.15p become 

5n(a l + a 2 ) , (l + T 2 n )5 n {Pl + Pl)\ (2nf/H n a r 6 n (2n) 2 b 2 n 2 , 



71 



+ 



(l + r n )3(l-r n ) 2(l + r n ) 5 (l-r, 



(l + r„)% 2 



2(1 + r n ) 



(5.16) 

Now we choose a n and b n so that the coefficients in the right hand side of (|5.16p tend to finite limits, 
say 1/2, giving 

(1 + Tnfol 



and 



(In 



2 7 /3 n 4/3^ ' 



V(l + Tn) 



2ny/S n 

Inserting into (I5.15p . we now choose 5 n so that we get a finite limit for the pre-factor, say 



(5.17) 
(5.18) 



2 1 / G ^(l + T n )n 2 / 3 a n b n e (^Hi-m) _ (1 + T „) 3 a[ 



n 3 / 2 a n 6 n 
It is easy to check that the choice 
s (l + r n ) 3 / 2 a n 



219/6^3/2 n5/3j 5/2 



exp 



2 1 / 3 n 4 / 3 <5 2 
(l + r„) 3 a\ 



2V6 V61ogcr n 



n 



2/3 ( 3 log a n - - log(6 log a n ) - log ( 2 1 / 4 vr(l + r, 



(5.19) 

V7T 



(5.20) 



satisfies (|5.19p . Note that this leads to a finite, non-zero limit for the kernel only on the diagonal, 
771 = r/ 2 ; if 771 ^ 772, we have + -^2)7(2^)1 > 1, implying that Mn^iQiA?) -> 0. By definition, 



;i+t„)« 1+ 



1 + r„ 



but since a n 3> 5 2 , we may expand in 5 n , neglect the second order term and choose 



(1 + O 1 + 



2(l+r n ) 2 / 

Substituting back the leading term of (|5.20p into (|5.17p and (|5.18p . we can thus choose a n and b n given 
by (12. lQh and (12. lip in the statement of the theorem. With these choices, we get that 

1 



the kernel of a two-dimensional Poisson process with intensity 7r 



-1/2^-^-5 



Again, it can be shown that there is a uniform integrable bound B^ ((^) = C^ e _ 2^+»? ) > 
{■M-n Tn ((, C)| for Re(C) > £o ; proving the convergence of the point process and last particle distri- 
bution to that of the limiting process. 
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6 Real integral representation of the interpolating Airy kernel 

6.1 Proof of Proposition 3] 

The proof of the alternative representation ()2.14j) is a simple matter of verification. Starting from 

(EH, 

— T'C^i+^i) f f r°° — 

iC^C^^)^-^— du dv dte it(«+.)-|<x 2 « 2 +|« 3 +iCi«-|<x 2 . 2 +|. 3 +if 2 ^ (61) 

4ir b / 2 a J r J r J 

where we have written the denominator as a convergent integral. Putting zj := Q + a 4 /A, the exponent 
in (16. ip can be trivially rewritten as 

/ \ ^~ 2 2 ^ S . j. ^22 ^ *-$ 

+ v) — a u H — u + iQ\u av H — v + iQ 2 v 

2 3 2 3 

, (7^ a 2 . 
=lt(u + V) ~ — + — [Zx + 2T 2 ) 

/ a 2 \ i ( a 2 \ 3 _ ( a 2 \ i ( a 2 ^ 3 

+ ZZi { u + 1 y) + 3{ u + 1 Y) + lZ2 \ v + l Y) + ?,\ v + l Y 

Using the well known integral representation of the Airy function, 

Ai(x) = — [ du e ixu+ ^ 

^ 2 2 

we obtain (j2. 14[) . upon substituting q = u + i 2 ^ and p = v + i 2 ^. 

6.2 Decay of density 

The asymptotic eigenvalue density when a n — > oo can be computed directly as the large n limit of the 
density, as expressed in (|2.18p , by an appropriate steepest descent analysis of /Cn T ™ • We will instead 
give a simple heuristic argument by calculating the large a limit of the density of the interpolating Airy 
process, after rescaling by a factor a, and recover the same expected decay of the density. Roughly 
speaking this amounts to saying that the same limit is obtained on the one hand if a n — > oo with n, 
and on the other if the large n limit is taken with a remaining constant, and then letting a tend to 
infinity. 

Changing variables t = av in the integral in (|2.14j) . we obtain 

2 — ct 6 +<t 3 £ — r}^ poo 

p a (0 := <r 2 lC*(a(, a() = a 66 / dv e - 3 M(Re ie )Ai(R e - ie ), 

V 71 " Jo 

where R = a\J (a 3 /A + v + £) 2 + r] 2 and tan# = r\(a 3 jA + v + £) . Using the large argument asymp- 
totics of the Airy function, (|2. 15|) . this can be written 

WO = ^Sr^r *' ^'- |R " lc °* (f »d + «(!))• 

Since 9 <C 1, we can expand in 6, and replace the highest order term with its constant value 9q = A-q/a 3 , 
leading to the approximate integrand e 9 "^"\ where 

g a {v) = va s - — 1 + 1 



6 V a 3 J \ a e 
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defined for v > — <t 3 /4 — £. For fixed £, g<j attains its maximum at i>o ~ — £, which stays of order 1 as 
a — > oo. Since g'l(v) is also of order 1 for w <C <r, and g'"{v) = 0(<r~ 3 ), expanding around v = up to 
second order in v will give the main contribution to the integral in (16. 2D . Neglecting terms of order 
a' 1 , we calculate ^(0) = -a 6 /6 - a 3 £ + rf - £ 2 , ^(0) = -2f , and ^(0) = -2, which gives 

„f°" 6 +o" 3 £ — r? roc k i roc i 

as <T — )• oo. 

7 Open problems 

There are several open problems left for future work. First of all, regarding universality, it would 
be very interesting to show that the interpolating kernel we have found also appears in the scaling 
limit for models with more general non-Gaussian weight functions, for instance the class of harmonic 
potentials. 

A second question is regarding the other symmetry classes f3 = 1,4, where we expect different 
interpolating Airy kernels, generalizing the corresponding Airy kernels for real eigenvalues. Perhaps 
the relations between the corresponding kernels pointed out recently in [5\ could open a simple way 
to obtain these. 

A somewhat more ambitious task is to investigate whether or not the relation between the Tracy- 
Widom distribution and the Painleve II solution has an analogue in the complex plane. It is possible 
that this will require a better understanding of the corresponding Riemann-Hilbert problem mentioned 
earlier. 
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